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Abstract 

Ground  Moving  Target  Indicator  (GMTI)  radar  maps  echo  data  to  range  and  range-rate, 
which  is  a  function  of  a  moving  target’s  velocity  and  its  position  within  the  antenna  beam 
footprint.  Even  stationary  clutter  will  exhibit  an  apparent  motion  spectrum  and  can 
interfere  with  moving  vehicle  detections.  Consequently  it  is  very  important  for  a  radar  to 
understand  how  stationary  clutter  maps  into  radar  measurements  of  range  and  velocity. 
This  mapping  depends  on  a  wide  variety  of  factors,  including  details  of  the  radar  motion, 
orientation,  and  the  3-D  topography  of  the  clutter. 
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Foreword 


Recently,  questions  have  arisen  over  the  precise  interrelationships  of  various  parameters 
and  characteristics  in  GMTI  radar  mode.  More  recently,  new  questions  have  arisen  over 
the  impact  of  target  scene  topography  on  GMTI  false  alarms.  The  need  to  account  for 
scene  topography  is  reasonably  well-known,  and  Sandia-designed  GMTI  systems 
routinely  include  provisions  to  do  so.  This  was  discussed  to  some  extent  in  a  section  of 
an  earlier  limited-release  report. 1 

This  report  details  the  analytical  basis  for  the  GMTI  range  and  velocity  measurements, 
and  how  clutter  is  perceived  by  the  GMTI  radar. 
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1  Introduction 


Ground  Moving  Target  Indicator  (GMTI)  radar  is  a  technique  whereby  a  radar  measures 
delay  and  delay-rate  from  a  sequence  of  pulses,  and  infers  a  moving  target’s  range  and 
velocity  from  these  measurements.  However,  the  transmogrification  of  the  radar 
measurements  into  useful  range  and  velocity  calculations  contains  subtleties  that  may 
result  in  errors  unless  these  subtleties  are  accommodated  and  the  otherwise  resulting 
errors  are  mitigated. 

This  requires  an  appreciation  for  the  finer  points  of  range  and  Doppler  interaction  for 
radar  echoes  from  moving  vehicles  in  various  radar  and  vehicle  geometries  and  motion. 

For  this  report,  we  will  limit  ourselves  to  an  exo-clutter  GMTI  mode,  that  is,  GMTI  from 
a  single  antenna  phase-center. 

In  the  subsequent  sections  we  address  a  number  of  issues,  including  (but  not  limited  to) 

•  the  accuracy  of  a  target  velocity  measurement, 

•  the  location  and  shape  of  clutter  in  the  range-Doppler  map, 

•  the  effects  of  antenna  rotation  on  target  scene  illumination,  and 

•  the  effects  of  target  area  topography  on  velocity  measurements. 

Examples  are  given  to  illustrate  these  issues. 

Clutter 


A  reasonable  question  is  “Why  is  the  location/shape  of  the  clutter  band  important?” 

First,  the  width  of  the  clutter  band  in  velocity-space  defines  the  Minimum  Detectable 
Velocity  (MDV)  for  an  exo-clutter  GMTI  radar.  Things  that  cause  the  clutter  band  to 
widen  or  narrow,  also  cause  the  MDV  to  change. 

Second,  the  clutter  band  denotes  a  ‘stay-out’  region  where  otherwise  clutter  is  expected 
to  generate  an  unacceptable  degree  of  false  alarm  detections.  As  the  clutter-band  moves, 
so  too  should  the  stay-out  region. 

Third,  the  location  of  the  center  of  the  antenna  beam  footprint  in  velocity-space  defines 
the  reference  velocity  for  any  moving  target,  that  is,  the  difference  between  this  velocity 
and  the  raw  measured  velocity  of  a  target  detection  represents  the  best  estimate  of  the 
target’s  own  range-rate.  As  the  center  of  the  beam  moves  in  velocity-space,  so  too  does 
the  reference  velocity. 
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" Nothing  happens  until  something  moves. 
—  Albert  Einstein 
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2  The  Range-Velocity  Map 


A  model  for  the  video  phase  history  data  for  a  stretch-processed2  GMTI  Coherent 
Processing  Interval  (CPI)  is  derived  in  Appendix  A.  It  is  expressed  as 


Xv(i,n)  =  Ar 


exp- 


-^-(t  -  Co  )fc,o,0  “  rc,0,0  fe 


2<yr,  o 


(v/os,s,0,0  vlos,c,0,o}^pn 


+ 


2  (c,0,0  rc,0,o)(v/c«,.?,0,0  vlos,c,0,o)Tpn 


(1) 


where 


{yios,s,  0,0  vlos,c,0,o)~ 


r.?(h)  +G,o,o) 

o' 

+ 

o 

v  r.?(h)  +G,o,o| 

rc(fi)  +G,0,o)|J 

•  V 


c  vlos,t,  0,0 


(2) 


and  individual  parameters  are  identified  as 

i  =  the  index  for  ADC  samples  of  a  single  pulse’s  echo  data, 
n  =  the  pulse  index  within  a  CPI, 

Ar  =  the  target’s  echo  strength. 

coT  o  =  the  transmitted  reference  radian  frequency, 

yT  0  =  the  transmitted  reference  radian  chirp  rate, 

=  a  factor  that  depends  on  line-of-sight  velocity  (see  Appendix  A), 
r?  (l  )  =  the  vector  difference  between  radar  location  and  target  location, 
rc(t)  =  the  vector  location  of  the  radar  with  respect  to  the  SRP, 
rs  n  J  =  the  nominal  range  between  radar  and  target, 
rc  n  j  =  the  nominal  range  between  radar  and  SRP, 
vc  =  the  radar  velocity  vector, 

vlos,s,n,i  =  the  nominal  line-of-sight  velocity  between  radar  and  target, 
vlos  cni  =  the  nominal  line-of-sight  velocity  between  radar  and  SRP, 
vlos  t  n  i  =  the  nominal  line-of-sight  velocity  component  of  just  the  target, 

Ts  =  the  ADC  sampling  period, 

Tp  =  the  pulse  period,  and 

c  =  the  velocity  of  propagation.  (3) 
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A  more  complete  description  of  these  parameters  can  be  found  in  Appendix  A. 

The  Scene  Reference  Point  (SRP)  is  the  location  at  which  the  antenna  is  pointed,  and  to 
which  any  motion  compensation  is  referenced.  It  defines  zero  range  offset,  and  zero 
Doppler  offset  at  that  range. 

The  third  phase  tenn  is  a  consequence  of  the  Residual  Video  Phase  Error  (RVPE)  term 
resulting  from  stretch  processing  (deramping  the  chirp  upon  receive).  It  is  typically 
ignored  for  GMTI  applications,  and  we  shall  do  so  hereafter. 

In  the  first  phase  term,  to  simplify  the  following  discussion,  we  define  the  range 
difference 

rsc  =ta,0,0-rc,0,o)-  (4) 


In  the  second  phase  tenn,  to  simplify  the  following  discussion,  we  define  the  line-of-sight 
velocity  difference 


vIos,sc  (v/os,s,0,0  Vos, c, 0,0  )  • 


(5) 


The  model  then  becomes 


Xv(i,n)=  Ar  exp] 


■— (1  -CoVsJsi 


2  CO' 


T,  0 


vlos,scTpn 


(6) 


An  examination  of  the  model  for  the  video  phase  history  data  shows  that  the  phase  is  a 
function  of  the  two  data  indices  n  and  i. 

A  Discrete  Fourier  Transforms  (DFT)  across  index  n  will  then  yield  an  Impulse  Response 
(IPR)  with  peak  at  apparent  velocity  offset 

v peak  ~  vlos,sc  (7) 


A  DFT  across  index  i  will  yield  an  IPR  with  peak  at  apparent  range  offset 

1  peak  ~  (l —  Co  Ysc  (8) 

We  now  examine  these  in  turn. 
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2.1  Moving  Target  Range 

2.1.1  Range  Difference  Expansion 

The  entity  of  concern  is  (l  )rSc  •  From  Appendix  A,  we  recall  that  the  factor  is 
given  by 

-  (v/os,s,0,0  -  vlos,c, 0,0 )  if  we  modulate  yT 

Co=jC  2  •  (9) 

-v/(W  o,0  if  we  don’t 

l  c 

For  all  reasonable  velocities  of  interest, 

vlos,s, 0,0  <<:  c  •  (10) 

Consequently,  while  perhaps  important  for  focusing  a  Synthetic  Aperture  Radar  (SAR) 
image,  we  may  presume  that  <^0  Fas  negligible  consequence  for  scaling  rsc  in  a  typical 
GMTI  mode.  That  is,  we  will  presume 

(l  ~~  Co  Ysc  ~  rsc  • 

We  identify  the  range  difference  as 

rsc  =  (g,0,0  ~rc,0,o)=  r.v(h)  +C,0,o|“  rc(?0  +C,0,o|  •  (H) 

We  make  use  of  the  following  equivalence 

+C,o,o)=rc(?o  +C,o,o)-so(h)  +C,o,o)-  (12) 

For  convenience  and  clarity,  we  recognize  that  these  are  constant  values  with  time,  and 
therefore  dispense  with  the  overt  time  arguments.  Consequently  we  now  write  this  as 

r5=rc-s0.  (13) 

We  note  that 

(14) 
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which  allows 


’sc 


f 

lii2  ) 

kd 

1  ,  lso|  2rc  •  s0 

i  i2  i  a 

U 

rT* 

(15) 


For  a  number  of  SAR  image  formation  algorithms  (e.g.  Polar  Format  processing3'4),  we 
use  the  Taylor  series  expansion 


-J 1  +  q  ~  1  + 


2  3 

a  a  a 

- 1 - h  ...  . 

2  8  16 


(16) 


and  then  often  only  the  linear  term.  Doing  so  allows  the  approximation 


'SC 


s0  rc*SQ 

2|rJ  IrJ 


(17) 


This  is  only  really  well  approximated  for  rsc  «  |rc| .  While  often  valid  for  SAR,  it  is  less 
tractable  for  the  large  swaths  of  a  wide-area  search  GMTI  mode. 

Nevertheless,  even  near  the  center  of  the  range  swath,  for  the  region  where  in  fact 
rsc  «  |rc| ,  we  will  observe  behavior  that  is  of  interest  to  us.  Consequently  we  will 

proceed  with  the  simplified  approximation  for  rsc  simply  to  gain  some  insight  to  the 
phenomenon. 

We  identify  the  target  coordinates  as 

So  =sxx  +  syy  +  szz  (18) 


where  x,y,z  identify  orthogonal  unit  vectors.  We  identify  the  radar  coordinates  as 

rc  =  —rc  cos y/0y  +  rc  sin  y/Qz  ( 1 9) 


where 


rc  =  nominal  scalar  range  between  radar  and  scene  center,  and 

y/Q  =  nominal  grazing  angle  at  scene  center.  (20) 

Consequently,  we  may  simplify 
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(21) 


rsc  *cos y/QSy  -sin i//qsz  +(yx2  +  sy2  +  sz2j/2rc  . 

At  this  point  we  make  some  observations. 

•  The  radar  measures  (l  -  )rsc  .  In  general,  rsc  depends  on  exclusively  target 
location.  There  is  a  minor  dependence  in  what  the  radar  measures  via  £g  ,  but  this 
is  a  very  minor  factor. 

•  In  general,  rsc  depends  on  offsets  from  the  SRP  in  all  three  dimensions,  namely 
all  components  sx,sy,sz.  However  the  dominant  dependence  is  on  the 
components  sy ,  sz  . 

•  Dependence  of  rsc  on  target  height  increases  as  the  grazing  angle  increases.  This 
is  the  familiar  range-layover,  or  foreshortening,  effect. 

To  first  order,  we  can  often  simplify  the  range  difference  model  even  more  to 

0  -  CoVsc  ~  -  r(|  *  |°  =  COSI//0.S'V  -  sin  y/Qsz  .  (22) 

The  fidelity  of  the  various  range  models  are  illustrated  in  Figure  1 . 

2.1.2  Range  Resolution 

The  resolution  with  which  we  can  measure  range  is  given  by  setting 

— (i-foAv.  =2* 

c  1=1 

rsc  ~Pr 

where 

pr  =  [slant]  range  resolution. 

This  is  solved  to  yield 

2  jvc  _  c 
Pr  “  (\-^)TsI2yTfi  =  (1-^)2 Beff 


(23) 


(24) 


(25) 
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where 


B. 


eff 


_  TsJyTfi 

In 


the  nominal  effective  bandwidth  in  Hz  contained  within  the 


data. 


(26) 


The  effective  bandwidth  Bejj  is  often  very  nearly  the  actual  transmitted  bandwidth,  but 
may  be  less,  especially  if  pulses  are  short  compared  to  the  range  swath  being  processed. 

Very  typically  £))  « 1 ,  so  it  is  often  ignored.  This  implies  that  an  adequate  expression 
for  range  resolution  is 


Pv 


(27) 


Note  that  this  is  the  achievable  resolution  of  the  system,  and  does  not  account  for  any 
coarsening  due  to  sidelobe  filtering  such  as  with  window  functions  prior  to  the  range 
compression  (the  DFT  across  index  i). 


200 
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2  100 

2 

i— 

CD 
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CD 

C 
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c 

03 


50 
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n  o 
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Figure  1.  Example  slant  range  error  of  various  models  for  rsc.  Swath  was  4  km  on  the  ground  at  10 
km  center  slant  range  from  4  km  altitude,  with  400  m  azimuth  offset. 
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2.2  Moving  Target  Velocity 

2.2.1  Line-of-Sight  Velocity  Expansion 


We  now  investigate  the  term 


vIos,sc 


>};los,s,  0,0 


vlos,c,  0,0 


^(h)  +^,o,o) 
v  rWo  +^cAo] 


rc(?Q  +tc,0  ,o) 

rc(t0  +tco,o'y 


•  V, 


vlos,t,  0,0  • 

(28) 


As  with  the  range  analysis,  we  make  use  of  the  following  equivalence 

rs(h)  +^,o,o)  =  rc(h)  +^,o,o)“so(h)  +  ^,o,o) 


(29) 


and  for  convenience  and  clarity,  write  this  as 
r?  =rc-s0. 


(30) 


Furthermore,  we  recall  that 


r  =  r  —  T 
'  SC  l1^  I  lc\- 


(31) 


This  allows  the  expansion  in  terms  of  just  the  radar  geometry  and  the  target  geometry  as 


vlos,sc 


f 

\ 

rc 

-1 

J 

rc  *  W 

s0*vc 

v  rc  rsc 

rc 

rc  +  rsc 

'  vlos,t, 0,0 


(32) 


which  can  be  rearranged  to  the  somewhat  more  convenient  form  for  our  purposes 


vlos,sc 


f  > 

1  1 

rc  *  vc 

1 

s0*vc 

H 

c°l- 

+ 

|rc 

1+Nc 

|rc| 

l  J 

V  rc  J 

'  vlos,t,0fi  ■ 


(33) 


Note  that  we  can  employ  the  Taylor  series  expansion 

— - —  =  1  -a +  o2 -a3 +o4 +...  =  V(-o),?  for  lal <  1 ,  and  n  >  0  .  (34) 

l  +  o 


-  15  - 


The  number  of  terms  required  for  high  fidelity  depends  on  just  how  much  a  <  1 . 
In  any  case,  this  allows  the  rearrangement  to 


vlos,sc 


I 


-1 


n  V  IV| J 


rc  •  vc 


s 


n  V  l*c| J 


S0*VC 


'  vlos,t, 0,0  • 


(35) 


The  relationship  of  rsc  to  Sq  was  explored  in  the  previous  section.  Recall  that  the  radar 


in  fact  measures  rsc  . 


While  for  wide-area  search  modes  we  may  often  have  large  range  swaths  compared  to  the 
center  range,  it  remains  instructive  to  examine  the  region  near  the  center  of  the  swath, 
where  rsc  <<  |rc  |  •  In  this  region,  the  model  for  line-of-sight  velocity  can  be  approximated 
as 


vlos,sc 


(  > 

r 

'sc 

rc  •  w 

so  *  A 

Uj 

rc 

rc 

vlos,t,0, 0  • 


This  can  be  expanded  to 


vlos,sc 


(  1  I2 

Pol 
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• 

1 

rc  •  vc 

s0*vc 

V^rc  2 

2 

rc  2 

kd 

|rc| 

vlos,t, 0,0  • 


(36) 


(37) 


At  this  point,  to  gain  further  insight,  we  need  to  pick  a  coordinate  frame,  and  constrain 
the  flight  path. 

We  shall  continue  with  target  coordinates  and  radar  coordinates  as  previously  defined. 
Accordingly,  we  identify  the  radar  velocity  vector  as 

W  =  ~va  sin  0S*  +  va  cos 0sy  +  W*  (38) 

where 

va  =  horizontal  (with  respect  to  target  plane)  radar  velocity  magnitude, 
vz  =  vertical  radar  velocity, 

0S  =  squint  angle  as  projected  into  target  plane.  (39) 


We  will  use  the  convention  for  squint  angle  as  depicted  in  Figure  2.  This  accounts  for 
the  negative  sign  in  front  of  the  x  component  in  the  radar  velocity  vector. 
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Putting  all  this  together  yields  the  expansion 


r  v a  sinffySy  Vacos^0cos^5x2  vz  sin  y/{)sx2  ^ 


vlos,sc 


2|r 


2|r 


2  2 
va  sin"  y/Q  cos 0ssy  va  cosy/0  cos 6ssy 


2  b 


vz  cos  y/Q  sin  y/(jsy  vz  sin  y/0sy 


bl  2|rc|2 

va  cos  y/Q  cos  6S  sin  y/QSz  va  cos  y/Q  cos  6ssz 


2  b 


vzszcosV0  vzsin^05z2 


2|r 


'  vlos,t,  0,0 


If  we  insist  on  straight  and  level  flight  ( v,  =0  ),  then  this  reduces  to 


(40) 
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f  2 

va  sin  Ossx  vqcos^Qcos  0ssx 


v I  os, sc 


2|r 


2  2 
va  sin  y/Q  cosGssy  va  cos  1//Qcos0ssy 


2  b 


Va  COS  ^0  COS  6S  sin  va  cos  V'O  cos  #Az 


2  b 


~vlos,t,  0,0 


(41) 


-2 


By  reducing  our  model  even  further,  by  neglecting  the  rc  tenns,  we  arrive  at 


f  .  2  n  h 

va  sin  6ssx  va  sin  i//0  cos  0ss} 


vlos,sc 


+ 


bl  bl 

va  cosi//q  cos0s  sin  i//0sz 
|rc- 1 

'  vlos,t,  0,0 


(42) 


At  this  point  we  make  some  observations. 

•  The  radar  measures  v/os%  s.c . 

•  In  general,  vios  sc  depends  on  both  the  target  velocity,  and  the  target  location. 

•  In  general,  vjos  sc  depends  on  offsets  from  the  SRP  in  all  three  dimensions, 
namely  all  components  sx,sy,sz. 

•  Dependence  of  vjos  sc  on  target  height  increases  as  the  squint  angle  approaches 

zero,  that  is,  as  the  radar  operates  towards  the  direction  of  the  velocity  vector  of 
the  radar. 

Remember  that  this  model  is  only  really  valid  near  the  center  of  the  range  swath. 
However,  we  do  expect  that  the  dependencies  on  geometry  will  not  lessen  appreciably  as 
distance  from  swath  center  increases. 
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Tarset  Height  Equivalent  Azimuth  Shift 


An  important  question  is  “How  sensitive  is  the  apparent  clutter  velocity  to  target  scene 
height?” 

We  answer  this  by  relating  line-of-sight  velocity  due  to  stationary  target  height,  to  line- 
of-sight  velocity  due  to  stationary  target  azimuth  offset.  That  is,  we  set 

va  sin  6ssx  _  va  cos^Q  cos#y  sin  y/Qsz  „ 

|rc|  lrc| 


Solving  this  yields 


sx  = 


^sin(2y/o)cot0s. 


(44) 


Note  that  this  is  dependent  on  both  squint  and  grazing  angles.  This  is  plotted  in  Figure  3. 
Note  that  at  small  angles,  the  relationship  can  be  much  more  than  one-to-one.  That  is,  a 
target  height  can  sometimes  be  equivalent  to  an  even  greater  target  azimuth  offset. 


0  10  20  30  40  50  60  70  80  90 

squint  angle  -  deg 


Figure  3.  Apparent  azimuth  shift  due  to  target  height  as  a  function  of  squint  angle,  for  various 
grazing  angles. 
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2.2.2  Velocity  Resolution 

The  resolution  with  which  we  can  measure  velocity  is  given  by  setting 


2(0TA) 

c 


VloS,Sc'Tpn 


\n=N 


=  2  TV 


v!os  ,SC'  =Pv 


where 


pv  =  velocity  resolution. 

This  is  solved  to  yield 

2  jvc  _  Aq 
Pv  ~  TpNlmrs,  ~2nTp 

where 

2  7ZC 

Aft  = -  =  the  nominal  transmitted  wavelength. 

®r,0 


(45) 


(46) 


(47) 


(48) 


We  note  that  NTp  is  the  total  collection  time,  that  is,  the  duration  of  the  CPI.  Note  also 

that  this  is  the  achievable  resolution  of  the  system,  and  does  not  account  for  any 
coarsening  due  to  sidelobe  filtering  such  as  with  window  functions  prior  to  the  velocity 
compression  (the  DFT  across  index  n ). 

Example 

Consider  a  GMTI  radar  that  operates  with  the  following  parameters. 

Aq  =  0.018  m  (Ku-band), 

N  =  256  pulses, 

Tp  =  250  ps  (PRF  =  4  kHz).  (49) 


These  parameters  allow  a  velocity  resolution  of 
pv  =  0.14  m/s. 


(50) 


Note  that  the  velocity  resolution  is  independent  of  radar  velocity.  We  do  note  that  in 
SAR,  sometimes  PRF  is  slaved  to  velocity,  making  velocity  resolution  dependent  on 
velocity  in  this  case.  However,  this  is  generally  not  true  for  GMTI. 
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2.3  Clutter  in  the  Range-Doppler  Map 

We  are  concerned  here  with  the  Doppler  signatures  of  stationary  clutter.  Of  course 
Doppler  yields  a  measure  of  relative  velocity.  So  we  limit  our  attention  to  vjossc  when 

vlos,t, 0,0  =  0  • 


Clutter  generally  provides  an  area  of  target  echo  responses.  The  only  limitation  is  the 
illumination  provided  by  the  antenna  beam  itself.  Consequently,  the  antenna  provides  to 
the  radar  echo  soundings  from  a  specific  region,  where  that  region  offers  a  spread  of 
velocity  indications. 

Clearly,  the  nature  of  the  spread  of  velocity  indications  is  defined  by  the  antenna  beam 
shape  and  orientation.  That  is,  Doppler  bandwidth,  and  hence  the  velocity  spread,  is 
limited  by  the  antenna  beamwidth,  primarily  the  azimuth  beamwidth. 

For  a  typical  elliptical  antenna  beam  pattern,  the  illuminated  region  on  the  ground  will 
also  be  elliptical  in  nature,  with  a  larger  region  illuminated  down-range  from  the 
boresight  than  from  ranges  closer  than  the  boresight.  Examples  will  be  given  in  later 
sections  of  this  report. 


We  begin  by  returning  to  the  exact  models  for  range  and  velocity,  namely 
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'  vlos,t, 0,0  ■ 


(51) 


Recall  that  vector  quantities  are  defined  as 


s0  =sxx  +  syy  +  szz, 

rc  =  -rc  cos  (//qj)  +  rc  sin  i//qz  ,  and 

vc  =  ~va  sin  0SX  +  va  cos 0sy  +  VZ^  •  (52) 

Consequently,  we  calculate  the  various  dot  products  as 
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(53) 


rc  •  s0  =  —rc  cos  i//0sy  +  rc  sin  y/QSz  , 

rc  •  vc  =  ~varc  cosi//()  cos  6^  +  vzrc  sin  i//() ,  and 

s0  •  vc  =  -va  sin  0ssx  +  va  cos  0ssy  +  vzsz . 


Since  GMTI  processing  nearly  always  involves  first  processing  a  range-Doppler  map, 
and  Doppler  is  proportional  to  line-of-sight  velocity,  we  are  interested  in  how  clutter 
maps  onto  a  grid  of  rsc  versus  viossc  . 

To  simplify  the  analysis,  and  to  make  some  specific  points,  we  will  confine  our  interest  to 
a  rectangular  scene  of  interest,  oriented  in  x  and  y  directions.  We  note  that  the 
illuminated  region  of  clutter  will  generally  not  be  a  rectangle,  as  previously  discussed. 
Nevertheless,  we  will  proceed  with  the  rectangular  scene  of  interest  and  observe  that  the 
effects  within  the  rectangular  scene  of  interest  will  apply  also  to  regions  outside  the 
rectangular  scene  of  interest. 

For  the  following  examples,  unless  noted  otherwise,  we  will  assume  the  following 
parameters 

va  =  100  m/s, 

|rc|  =  10  km, 

y/ o=30  degrees.  (54) 


2.3.1  Flat  Terrain  in  Straight  and  Level  Flight 

We  first  consider  a  flat  target  plane  with  sz  =  0 .  Straight  and  level  flight  implies  vz  =  0  . 


Example 

Consider  a  target  scene  with  stationary  reflectors  arranged  with  ground  truth  as  in  Figure 
4.  Now  consider  the  scene  being  interrogated  from  a  number  of  different  squint  angles. 
The  resulting  range-velocity  maps  are  depicted  for  squint  angles  of -90,  -60,  and  -30 
degrees  respectively  in  Figure  5,  Figure  6,  and  Figure  7. 

We  observe  the  following. 


•  Except  at  broadside  squint  angles,  there  is  a  range-dependent  shift  in  apparent 
clutter  velocity.  It  is  non-linear  with  range,  and  manifests  as  a  ‘comma’  shape. 

•  The  scene’s  velocity  band  narrows  in  velocity  as  squint  angle  moves  away  from 
broadside. 
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Figure  4.  Example  ground  truth  for  stationary  reflectors. 
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Figure  5.  Range-  velocity  map  for  -90  degree  squint  angle. 
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Figure  6.  Range-  velocity  map  for  -60  degree  squint  angle. 


20 


2500 

2000 

1500 

1000 

500 
E 

■  0 

o 

C/3 

-500 
-1000 
-1500 
-2000 
-2500 

-20  -15  -10  -5  0  5  10  15  20 


Vsc  -  m/s 


Figure  7.  Range-  velocity  map  for  -30  degree  squint  angle. 
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2.3.2  Elevated  Terrain  in  Straight  and  Level  Flight 


We  now  concern  ourselves  with  the  condition  sz  ^  0  . 

Example 

With  respect  to  the  previous  example,  we  now  raise  the  target  plane  by  200  m. 

Otherwise  it  remains  flat.  Now  consider  the  scene  being  interrogated  from  the  same 
squint  angles  as  those  of  the  previous  examples.  The  resulting  range-velocity  maps  are 
depicted  for  squint  angles  of -90,  -60,  and  -30  degrees  respectively  in  Figure  8,  Figure 
9,  and  Figure  10. 

This  exhibited  shifting  is  ‘layover’  in  SAR  images,  and  happens  in  both  range  and 
velocity,  depending  on  geometry.  As  a  sanity  check,  we  note  that  on  the  left  side  of  the 
flight  path,  and  forward  of  broadside,  as  the  stationary  target  is  elevated,  it  achieves  a 
smaller  angle  with  respect  to  the  velocity  vector,  hence  exhibits  a  more  positive  closing- 
rate,  and  hence  a  more  negative  range-rate,  that  is,  a  more  negative  line-of-sight  velocity. 
It  also  exhibits  a  closer  range.  Hey,  it  checks  out. . . 


Figure  8.  Range-  velocity  map  for  -90  degree  squint  angle. 
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Figure  9.  Range-  velocity  map  for  -60  degree  squint  angle. 
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Figure  10.  Range-  velocity  map  for  -30  degree  squint  angle. 
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2.3.3  3-D  Topography  Effects  -  Range  Dependent  Topography 

We  now  concern  ourselves  with  the  condition  sz  ^  0  ,  where  elevation  varies  with  range. 
The  radar  is  still  presumed  to  fly  straight  and  level. 

Example 

Now  consider  a  target  scene  with  stationary  reflectors  arranged  with  ground  truth  that 
includes  height  topography  as  in  Figure  11.  Now  consider  the  scene  being  interrogated 
from  a  number  of  different  squint  angles  matching  those  of  the  previous  examples.  The 
resulting  range-velocity  maps  are  depicted  for  squint  angles  of -90,  -60,  and  -30  degrees 
respectively  in  Figure  12,  Figure  13,  and  Figure  14. 

We  observe  the  following. 

•  Variations  in  target  height  will  cause  variations  in  apparent  target  velocity. 

•  An  elevation  offset  that  is  constant  at  any  one  range  will  cause  a  shift  in  the 
scene’s  velocity  band. 

•  The  degree  of  target-height  induced  velocity  variations  depends  on  the  squint 
angle,  increasing  as  the  viewing  geometry  moves  away  from  broadside. 


2.3.4  3-D  Topography  Effects  -  Azimuth  Dependent  Topography 

In  the  previous  section,  we  concerned  ourselves  with  a  range-dependent  clutter  elevation 
variation.  We  now  concern  ourselves  with  an  azimuth-dependent  elevation  variation  of 
the  clutter.  The  radar  is  still  presumed  to  fly  straight  and  level. 

Example 

Consider  the  same  GMTI  radar  as  in  the  previous  examples.  Now  consider  a  target  scene 
with  stationary  reflectors  arranged  with  ground  truth  that  includes  height  topography  as 
in  Figure  15.  Figure  16  illustrates  the  range-Doppler  map  when  interrogated  at  a  -30 
degree  squint  angle. 

We  observe  that  the  scene’s  velocity  band  has  widened. 

With  other  slopes,  the  scene’s  velocity  band  may  in  fact  narrow. 
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Figure  11.  Example  ground  truth  for  stationary  reflectors,  but  with  height  variations  of  400  m  peak- 
to-peak. 


Figure  12.  Range-  velocity  map  for  -90  degree  squint  angle. 
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Figure  13.  Range-  velocity  map  for  -60  degree  squint  angle. 
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Figure  14.  Range-  velocity  map  for  -30  degree  squint  angle. 
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Figure  15.  Example  ground  truth  for  azimuth  dependent  elevation  variation. 
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Figure  16.  Range-  velocity  map  for  -30  degree  squint  angle. 
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2.3.5  Flat  Terrain  -  Radar  Vertical  Velocity 


In  the  previous  sections,  the  radar  was  presumed  to  fly  straight  and  level.  Here  we  depart 
from  this  presumption,  and  now  presume  a  vertical  velocity  component. 

Example 

Consider  the  same  GMTI  radar  as  in  the  previous  examples,  except  that  the  radar  now 
has  a  vertical  velocity.  Reconsider  the  flat  target  scene  with  stationary  reflectors  arranged 
with  ground  truth  as  in  Figure  4.  Figure  17  illustrates  the  range-Doppler  map  when 
interrogated  at  a  -30  degree  squint  angle. 

We  observe  that  the  scene’s  velocity  band  has  shifted  as  a  function  of  range,  just  due  to 
the  radar  vertical  velocity. 


Figure  17.  Range-velocity  map  for  -30  degree  squint  angle. 


-31  - 


2.4  Antenna  Beam  Alignment 


The  position  and  orientation  of  the  clutter  region  in  a  range-velocity  map  will  also 
depend  on  just  what  clutter  is  illuminated  by  the  radar.  For  example,  if  the  antenna  were 
to  exhibit  an  unknown  pointing  bias  in  the  direction  of  the  radar  velocity  vector,  then  the 
clutter  band  would  be  shifted  in  the  range-velocity  map  towards  more  negative  velocities 
(range-rates).  Remember  that  the  most  negative  range-rate  is  in  the  direction  of  the  radar 
velocity  vector,  as  this  is  the  direction  that  range  is  decreasing  the  fastest. 

Most  radar  systems  have  antennas  that  exhibit  elongated  antenna  beam  patterns,  often  a 
fan  beam,  where  elevation  beamwidth  is  greater  (sometimes  substantially  so)  than  the 
azimuth  beamwidth.  A  three-axis  gimbal  (or  equivalent)  will  typically  allow  adequate 
beam  rotations  to  a  proper  desired  orientation  and  direction,  most  often  with  maximal 
symmetry  in  pan  and  tilt  directions. 

If  using  a  two-axis  gimbal,  often  the  missing  degree  of  freedom  is  one  that  allows  for 
rotation  around  the  boresight  of  the  antenna  beam,  sometimes  called  “twist”.  Clearly,  a 
rotation  around  the  boresight  for  a  fan-beam  will  cause  a  corresponding  rotation  of  the 
footprint  of  the  antenna  beam  on  the  ground. 

This  of  course  defines  the  clutter  that  is  illuminated,  and  hence  the  basic  nature  of  how 
clutter  appears  in  a  range-velocity  map,  quite  independent  of  and  in  addition  to  any  other 
factors  that  affect  the  clutter  map. 

We  now  investigate  the  effects  of  a  rotation  of  an  elliptical  antenna  beam. 

2.4.1  The  Antenna  Beam  Pattern  Model 


Consider  an  antenna  beam  where  we  define  the  angular  ‘edge’  as  an  ellipse.  This  is 
illustrated  in  Figure  18.  The  equation  that  governs  this  shape  is 


<t>h  !  A2  =1 

[daz/lf  (%/2)2 


(55) 


where 


0az  =  nominal  azimuth  beamwidth  of  the  antenna  beam, 

0el  =  nominal  elevation  beamwidth  of  the  antenna  beam, 

(f>l%  =  horizontal  offset  of  beam  edge,  and 

<f>v  =  vertical  offset  of  beam  edge.  (56) 
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With  some  foresight,  we  also  identify  derivative  parameters 


tan  6  =  — 

<K 


and 

p2=42+42. 


(57) 


(58) 


In  these  terms,  we  can  modify  the  ellipse  equation  describing  the  antenna  beam  pattern  to 


p1  sin2  0  +  p2  cos2  6  _  ^ 

(ejif  +  fe/2  )2 


Figure  18.  Model  for  antenna  beam  shape. 
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We  may  conveniently  rotate  the  ellipse  by  introducing 
a  =  antenna  beam  rotation  about  its  boresight. 


(60) 


This  is  illustrated  in  Figure  19.  This  causes  a  modification  of  the  antenna  pattern  ellipse 
equation  to 


p1  sin2 (d-  a)  p1  cos2(#-a)  _ 

(<W2)2  +  fc/2  f 


(61) 


However,  it  remains  true  that 
(f>h  =  p  sin  6 

<j)v  =  p  cos  6 .  (62) 


Figure  19.  Model  for  rotated  antenna  beam  shape. 
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2.4.2  The  Target  Scene  Projection  Model 


We  will  assume  that  the  antenna  boresight  is  directed  towards  the  center  of  the  scene 
being  interrogated.  The  question  then  becomes  “How  do  angular  offsets  from  the 
antenna  boresight  project  onto  the  target  plane?”  We  will  further  assume  a  horizontal  flat 
target  plane  with  no  elevation  offset.  This  is  illustrated  in  Figure  20. 


radar 


It  is  straightforward  to  derive 
..  _  rc  sin  i//q  tan 4 

sin(^o-^v) 


(63) 


and 


y  =  rc  sin  y/Q  (cot(^0  -  &  )  -  cot  y/Q ) . 


(64) 


These,  with  the  antenna  beam  shape  model,  will  allow  us  to  plot  projections  of  the  beam 
shape  onto  the  ground.  We  do  so  next. 
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2.4.3  Target  Scene  Projection  Examples 


What  follows  are  several  examples  of  projecting  rotated  elliptical  antenna  beams  onto  a 
horizontal  flat  ground.  In  all  cases  we  will  presume 

6az  =  3  degrees, 

0el  =7  degrees, 
rc  =  10  km,  and 

y/ o=10  degrees  (65) 


We  illustrate  antenna  rotations  about  its  boresight  of  0  degrees,  30  degrees,  60  degrees, 
and  90  degrees  respectively  in  Figure  2 1  through  Figure  24.  We  make  several 
observations. 


•  Even  for  a  perfectly  oriented  antenna,  the  ground  projection  is  asymmetrical  in 
the  down-range  (y)  direction. 

•  A  rotation  of  the  antenna  beam  results  in  a  rotation  of  the  beam  footprint  on  the 
ground,  but  the  amount  of  rotation  of  the  projection  is  not  identical. 

•  For  a  vertical  fan  beam,  the  rotation  for  small  angles  is  in  the  same  direction. 

•  For  a  vertical  fan  beam,  antenna  rotation  also  results  in  a  smaller  range  swath 
being  illuminated. 

•  Even  a  horizontal  fan  beam  may  project  longer  in  the  down-range  direction  than 
in  the  cross-range  (azimuth)  direction.  This  is  because  of  the  grazing  angle. 

Figure  25  illustrates  the  relationship  of  rotation  in  the  ground  plane  projection  to  the 
rotation  of  the  antenna  beam  about  its  boresight  for  several  grazing  angles.  Note  that 
sensitivity  for  small  grazing  angles  decreases  with  grazing  angle.  While  an  interesting 
and  important  characteristic,  it  is  deceptive  to  monitor  only  this  and  ignore  how  the 
dimensions,  particularly  the  aspect  ratio,  of  the  projection  also  changes  with  antenna 
rotation. 
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Figure  21.  Antenna  pattern  and  ground  projection  for  0  degrees  of  antenna  rotation. 
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Figure  22.  Antenna  pattern  and  ground  projection  for  30  degrees  of  antenna  rotation. 
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Figure  23.  Antenna  pattern  and  ground  projection  for  60  degrees  of  antenna  rotation. 
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Figure  24.  Antenna  pattern  and  ground  projection  for  90  degrees  of  antenna  rotation. 
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Figure  25.  Ground  projection  rotation  versus  antenna  rotation  at  various  grazing  angles. 
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2.5  Clutter  Motion 


Up  until  now,  we  have  assumed  that  those  target  scene  features  that  we  normally  assign 
the  label  as  “clutter”  were  stationary  during  the  course  of  the  data  collection.  This  might 
be  stretching  the  truth  somewhat.  In  fact,  the  clutter  itself  may,  actually  does,  have  a 
motion  spectrum  primarily  due  to  wind.  Quite  simply,  wind  makes  clutter  move.  This 
motion  will  broaden  the  Doppler  spectrum  of  the  clutter,  and  broaden  the  velocity 
spectrum  that  clutter  occupies  in  a  range-velocity  map. 

Skolnik5  compiles  results  of  other  authors  to  indicate  that  windblown  clutter  exhibits  a 
Gaussian  distribution  that  is  characterized  by  its  velocity  standard  deviation,  <yv .  The 
table  that  Skolnik  presents  includes 

<jv  up  to  0.32  m/s  for  wooded  hills  with  40  kt  wind, 
crv  up  to  1.1  m/s  for  sea  echoes, 
crv  up  to  1 .2  m/s  for  chaff  in  25  kt  wind,  and 
crv  up  to  4.0  m/s  for  rain  clouds. 

Note  that  1  kt  ~  0.5 14  m/s. 

The  discussion  omitted  mention  of  mean  velocity  for  chaff  and/or  rain.  Things  that  move 
(e.g.  chaff,  rain,  trash,  etc.)  will  probably  have  a  significant  mean  horizontal  velocity  that 
may  be  in  the  region  of  the  wind  mean  velocity.  Presumably,  the  mean  velocity  for 
wooded  hills  was  near  zero,  since  trees  don’t  typically  translate  or  relocate  (except  for 
very  big  winds).  In  any  case,  we  would  expect  that  foliage  velocity  is  generally  limited 
to  wind  velocity,  and  perhaps  more  strictly,  foliage  velocity  is  likely  limited  to  wind 
velocity  variations  (e.g.  gusts,  etc.).  This  is  likely  also  the  case  with  cultural  clutter  items 
such  as  hanging  objects.  Windmills,  fans,  and  turbines  are  a  different  story. 

Nathanson6  shows  data  that  suggests  that  the  -3  dB  velocity  bandwidth  of  sea  clutter  is 
related  to  Hydrographic  Sea  State,  S,  by  the  approximation 

5V  *!  m/s.  (66) 
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2.6  Summary  Discussion 


The  location,  width,  and  shape  of  the  clutter  in  a  range-velocity  map  are  influenced  by 
any  of  a  number  of  scene  and  radar  flight  parameters,  including 

Radar  horizontal  velocity 
Radar  squint  angle 
Radar  vertical  velocity 
Radar  antenna  beam  footprint 

Target  scene  elevation  offset 
Target  scene  elevation  variations  in  range 
Target  scene  elevation  variations  in  azimuth 
Target  scene  composition  and  wind  effects 

An  accurate  prediction  of  how  the  clutter  will  map  into  range-velocity  space  will  depend 
on  all  of  these. 


2.7  Clutter  Mitigation 

In  exo-clutter  GMTI,  clutter  is  a  source  of  false  alarms.  To  mitigate  the  effects  of  clutter 
on  false  alanns,  the  GMTI  needs  to  ignore  detections  in  range-velocity  regions  that 
correspond  to  clutter.  For  this  to  occur,  the  radar  needs  to  know  (or  calculate)  where 
clutter  topography  maps  into  range-velocity  space. 

Three  general  classes  of  techniques  may  be  employed  to  mitigate  topography  effects. 

1 .  The  radar  is  given  topography  infonnation  to  estimate  the  location  of  clutter  in 
range-velocity  space.  This  may  be  in  the  fonn  of  Digital  Terrain  Elevation  Data 
(DTED)  input  from  a  library.  The  radar  also  needs  antenna  attitude  information. 

2.  The  radar  makes  its  own  topography  measurements,  perhaps  in  a  manner  similar 
to  Interferometric  Synthetic  Aperture  Radar  (InSAR  or  IFSAR).  This  requires  a 
second  antenna  phase  center  in  elevation.  The  radar  still  also  needs  antenna 
attitude  information. 

3.  The  radar  estimates  whether  a  pixel  in  the  range-velocity  map  belongs  to  clutter 
or  not  based  on  the  range-velocity  map  itself.  This  is  a  data-driven  adaptive 
technique.  Basically,  the  radar  tries  to  ‘detect’  the  clutter  and  identify  it  as  such. 

In  all  cases,  accurate  and  precise  radar  motion  information  is  required. 

Of  course,  combinations  and  variations  of  these  basic  ideas  may  also  be  useful. 
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“Get  your  facts  first,  then  you  can  distort  them  as  you  please.  ” 

—  Mark  Twain 
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3  Conclusions 

The  following  points  are  worth  repeating. 

•  In  exo-clutter  GMTI  radars,  the  ground  clutter  generates  false  alarms.  This  is 
undesirable. 

•  GMTI  radars  need  to  ignore  detections  in  regions  of  clutter.  To  do  this,  the  radar 
needs  to  know  how  clutter  maps  into  range-velocity  space  for  the  radar. 

•  The  nature  of  how  clutter  maps  into  range-velocity  space  for  GMTI  depends  on  a 
number  of  factors,  including  generally  radar  motion  and  viewing  geometry,  and 
the  topography  of  the  target  scene. 

•  The  radar  can  determine  the  clutter  characteristics  in  its  range-velocity  space  by 
several  different  techniques.  These  include  using  radar  motion/attitude 
information  with  target  scene  topography,  or  perhaps  by  trying  to  ‘detect’  the 
clutter  based  on  the  data  itself.  Additionally,  the  radar  may  measure  the  target 
scene  topography  using  other  techniques  such  as  IFSAR.  Combinations  of  these 
may  also  be  used. 
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Figure  26.  Christian  Andreas  Doppler  (29  November  1803  -  17  March  1853)  was  an  Austrian 
mathematician  and  physicist.  The  Doppler  effect  is  named  after  him.  The  Doppler  effect  is  an 
apparent  change  in  frequency  and  wavelength  due  to  relative  motion  between  a  wave’s  source  and 
an  observer. 
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Appendix  A  -  The  Phase  History  Model 


In  this  appendix  we  develop  a  model  for  the  radar’s  phase  history  data  as  collected  by  a 
GMTI  radar  using  stretch  processing. 

We  will  also  presume  that  the  GMTI  radar  dwells  on  a  reference  location  for  the  duration 
of  its  Coherent  Processing  Interval  (CPI).  We  name  this  reference  location  as  the  Scene 
Reference  Point  (SRP)  and  designate  it  as  the  origin  of  the  coordinate  frame  used  for  data 
collection  and  processing. 

We  apologize  up  front  for  this  being  an  exercise  in  keeping  track  of  all  the  various 
subscripts. 

The  Radar  Signals 

In  general,  the  radar  will  collect  data  over  a  CPI  using  N  pulses.  We  allow  the  pulse 
index  n  to  vary  as 


-N 

2 


^  N 
<n<  — . 

2 


(Al) 


At  some  range  from  the  SRP  the  radar  emits  a  pulse  towards  the  SRP  which  we  model  as 


XT{t,n)  =  rect 


expty 


f  Y 

<h,n  +(°T,n(t~tn)+^L{t~tn)2 


(A2) 


where 


tn  =  the  reference  time  for  pulse  n  , 

T  =  the  transmitted  pulse  width, 

(pT ,n  =  reference  phase  for  transmitted  pulse  n, 

coT  n  =  the  reference  radian  frequency  for  transmitted  pulse  n, 

Yt  n  =  the  radian  chirp  rate  for  transmitted  pulse  n.  (A3) 


We  define  the  pulse  envelope  function 


rec 


t(xHo 


|x|  <1/2 

else 


(A4) 


The  received  echo  signal  is  a  delayed  and  attenuated  version  of  the  transmitted  signal. 
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XR(t,n)=  Ar  rect 


(A5) 


where 

ts  =  the  echo  delay  time  for  pulse  n, 

Ar  =  the  modified  signal  amplitude  of  the  received  echo.  (A6) 


Recall  that  we  are  generally  presuming  a  relative  motion  between  radar  and  target, 
consequently  any  echo  delay  time  will  itself  be  time-dependent.  While  not  explicitly 
shown,  ts  is  itself  a  function  of  time. 


We  will  presume  that  our  GMTI  uses  stretch  processing,  whereby  the  received  echoes  are 
‘dechirped’  with  a  Local  Oscillator  signal  modeled  as 


XLO(t,n)  =  rect 


( t  t n  hn ) 


exPi  J 


J 


$L,n  +  ®L,tS$  hi  hw )  2  hi 


(A7) 


where 


tm  =  the  Local  Oscillator  reference  offset  time, 

Tl  =  the  Local  Oscillator  pulse  width, 

</>Ln  =  the  reference  phase  for  Local  Oscillator  pulse  n, 

coL  n  =  the  reference  radian  frequency  for  Local  Oscillator  pulse  n, 

yL  n  =  the  radian  chirp  rate  for  Local  Oscillator  pulse  n.  (A8) 


Application  of  the  Local  Oscillator  signal  to  the  received  echo  yields  a  video  difference 
signal  of  the  form 
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Xv(t,n)  =  Ar  rect 
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+  C0T,n(t~tn  ~ts)~C0L,n(t~tn  ~tm) 

+  ^L{t-tn  ~tsf  -^L(t-tn  -fmf 


(A9) 


This  signal  is  sampled  by  the  Analog  to  Digital  Converter  (ADC)  at  times 


t  -tn+  tm  +  Tsi 


(A10) 


where 


Ts  =  the  ADC  sample  spacing 
i  =  the  ADC  sampling  index  within  a  single  pulse. 

We  allow  for  a  multitude  of  samples  totaling  /,  with  index  i  such  that 


(All) 


(A  12) 


At  this  point  we  observe  that  ADC  sample  times  span  an  interval  of  TSI .  As  previously 

noted,  during  this  time  both  target  and  aircraft  are  moving.  Consequently  the  target  delay 
is  in  fact  changing  during  the  time  we  are  sampling,  as  well  as  on  a  pulse  to  pulse  basis. 
We  now  explicitly  show  ts  as  a  function  of  indices  i  and  n,  and  model  this  as 
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s,n,i 
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(A  13) 


The  resulting  video  signal  is  then  described  by 
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We  manipulate  this  to 


&T,n  ~  ^L,n 
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As  a  part  of  real-time  motion  compensation,  we  equate 
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(A  15) 


(A  16) 


This  simplifies  the  video  signal  model  to 


Xv(i,n)=Ar  exp]  j 
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(A  17) 


In  SAR  we  often  modulate  frequency  and  chirp  rate  to  facilitate  mitigation  of  migration 
tenns,  which  is  necessary  for  fine  resolution  imaging.  For  GMTI,  CPIs  are  short  enough 
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and  range  resolution  is  typically  coarse  enough  that  migration  is  not  a  significant 
problem. 


In  a  manner  similar  to  target  echo  delay  time,  the  echo  delay  time  to  the  SRP  is  also  time 
dependent  and  can  be  expressed  as 


t c  t c,n,i  t c,n,0 
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d(TJ ) 


c,n,i 
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(A  18) 


We  note  in  the  first  line  of  the  phase  tenn  the  overt  dependence  on  the  time-dependent 
echo  delay  time,  namely  what  we  will  identify  later  as  a  dependence  on  the  line-of-sight 
velocity.  We  might  to  advantage  compensate  this  term  somewhat  by  modulating 
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d(T,i)‘c- 


(A19) 


for  some  constant  yT  0 .  This  would  yield  a  video  signal  of  the  fonn 
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or  approximately 
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with  a  substantially  reduced  influence  by  the  radar  velocity.  Nevertheless,  we  will 
presume  constant  reference  frequency  and  allow  for  a  chirp  rate  modulated  in  this  fashion 
or  not  by  defining  a  parameter  such  that 
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(A22) 


Sandia-designed  GMTI  modes  designed  to  date  do  in  fact  modulate  the  chirp  rate  in  this 
fashion.  The  model  then  becomes 


-  \(0T  Q  +  yT, o(l  Cn  hi  h 

Xv(i,n)=  Ar  exp- 

j 

v  2  J) 

(A23) 


The  term  in  the  second  line  of  the  phase  of  this  expression  is  the  Residual  Video  Phase 
Error  (RVPE)  tenn.  It  can  frequently  be  mitigated  with  processing,  but  is  often  ignored 
for  coarse  resolution  processing,  too.  We  will  carry  it  along  for  now. 

Radar  Geometry  and  Motion 

In  general,  the  target  location  and  the  radar  location  will  both  vary  with  time,  we  define 
these  time-dependent  locations  as 

s(t)  =  the  vector  describing  the  location  of  the  target,  and 

rc\t)  =  the  vector  describing  the  location  of  the  radar.  (A24) 

These  are  illustrated  in  Figure  27. 
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Time-varying  quantities  will  vary  with  both  indices  n  and  i.  We  define  time -varying 
vector  quantities  that  depend  on  these  indices  as  follows. 

sni  =  the  vector  describing  the  location  of  the  target  for  pulse  n  and  sample  i, 
rc  n  j  =  the  vector  describing  the  location  of  the  radar  for  pulse  n  and  sample  i. 

(A25) 


We  identify  velocities  as  the  respective  time  derivatives  of  location,  namely 
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d_ 

dt 

d_ 

dt 


s(t)  =  target  velocity  vector,  and 
rc(t)  =  radar  velocity  vector. 


(A26) 


We  will  assume  these  velocities  to  be  constant  during  a  CPI. 

Furthermore,  we  identify 

r?(t)  =  rc(t)-s(t)  =  vector  of  radar  location  relative  to  target, 

v,  =  vc  -  V;  =  radar  velocity  relative  to  target.  (A27) 

Consequently,  we  can  equate 

r,  M  =  r,  ifn  +  tm  +  Tsi)+  (t  -  tn  -  tm  -  Tsi )ys  , 

=  rc(t„  +  tm  +  Tsi)  +  (t  -  tn  -  tm  -  Tsi)\c ,  and 

s»W=siife  +tm  +Tsi)+(t-tn -tm-Tsifrf  (A28) 


We  define  range  as  a  scalar  equal  to  the  magnitude  of  location  vectors.  That  is 
fvM  =  lrv(c  +  tm  +  Td]  =  range  to  the  target,  and 

rc  n.i  =  | rc(t„  +  tm  +  TJ)\  =  range  to  scene  center.  (A29) 


The  following  scalar  line-of-sight  velocities  are  also  useful 
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Over  short  CPIs  and  for  expected  target  vehicle  motions,  we  will  presume  these  to  be 
adequately  modeled  as  constants.  In  general,  vios  c(t)  is  known,  but  vlos  s (l )  and  vlos  t  (t ) 

are  not. 

We  note  that 

vfo^(0  =  -l44#(vc-vJ.  (A31) 

r,(0| 


This  relationship  can  also  be  rewritten  and  expressed  as 
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At  particular  sample  times, 
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Radar  Echo  Delays 

Echo  time  delays  are  related  to  range  by  the  velocity  of  propagation.  Specifically, 

ts  =  —  |rv(/  )| ,  and 
c 

4=-MT  (A35) 
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where 


c  =  the  velocity  of  propagation. 


(A3  6) 
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In  terms  of  sample  times, 


h,n,i  =  ~rs,n,i >  and 
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rc,n,i  —  vlos,c,n,i  ■ 
C 


It  is  adequate  to  model 


vlos,s,n,i  ~  vlos,s,n,0  ’  a^d 
vlos,c,n,i  ~  vlos,c,n,0  • 


(A37) 


(A3  8) 
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Our  model  for  the  video  phase  history  data  remains 


(  +  ^r,o(l  bbz  ^dKfyn,  0  bn) 

Xv(i,n)=  Ar  exp< 

j 

[+r2°fc.«.0  bn)2  jj 

where  we  now  identify 


(yios,s,n,  0  vlos,c,n,  0  ) 

C 


vlos,s,n,0 


if  we  modulate  „ 
if  we  don't 


(A40) 
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to  Pulse  Variations 


We  now  turn  our  attention  to  the  target  delay  characteristics  and  how  they  change  from 
pulse  to  pulse.  The  target  delay  tenn  ts0  n  is  modeled  to  be  static  during  a  pulse,  but 

varies  from  one  pulse  to  the  next.  By  assuming  a  constant  Pulse  Repetition  Frequency 
(PRF)  we  can  model  this  as 


2 

c 


t s,n,0  rs,n,  0 


(fy0,0  VloSySfifi'Tpft) 


(A42) 
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where 

rs  ().o  =  nominal  target  range  at  center  of  CPI,  and 
Tp  =  the  pulse  period  =  1/PRF. 

We  furthermore  identify  the  SRP  delay  as 

C,«,0  ~~rc,n,0  (rc,0,0  vlos,c,0,0^1pn) 

C  C 

where 

rc  o  o  =  nominal  range  to  the  SRP  at  center  of  CPI. 


(A43) 
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(A45) 


We  choose  to  operate  the  radar  so  that  the  Local  Oscillator  reference  time  tracks  the 
delay  from  the  scene  center,  namely 

=  tc,n, o  •  (A46) 


We  consequently  calculate  the  relative  target  delay  as 

(fs',«,0  —  b??)  —  ~rc,n, o)—  (f?, 0,0  —  ^c,0,o)^  i};los,s, 0,0  —  vlos,c,0,o)^pn  • 


(A47) 


Combining  this  into  the  video  phase  history  model  yields 


r  r 


Xv(i,n)=  Ar  exp 


\f°T  ,o  +  Yt.q  (l  ~  Cn  n 


\rs, 0,0  ~rc, 0,0  ) 


,  ivlos,s,(),()  vlos,c,0fi)^pn 

V  C 


^^i~{rs,0,0~rc,0,0)+~{vlos,s,0,0~vlos,c,0,0yPn 


Y2 


(A48) 


We  parse  the  various  phase  terms  to  yield 
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Xv (/, n)  =  Ar  exp 


— ®r,o(^,o,o  _rc,o,o) 
c 

2 


777)0  “  Cn  ){rs,0,0  “  rc,0,oKz 
®7\0 (l7o,y,,y,0,0  “  vlos,c,0,o)^pn 
7r,o0  ^ nfxios,s,o,o  ~  vios,c,o,o)r piiTsi 


~  (  Oy,0,0  ^c,0,o)~*"  (y/o,y,,y,0,0  vlos,c,0,o)^p^ 

2  1  c 


7  7 


or 


Xv(i,n)=Ar  exp] 


2<» 


'r,o  /  \ ,  2/7,0  /  \2 

■ — 0,0  -  rc,0,0  j+ 2  Vs, 0,0  -  rc,0,of 


+ 


+ 


2 
V  c 


-^■O  -  Cz  fe,0,0  -  rc,0,0  fei 

c 

2(Ot,o  i  7., 

ly/os,.y,0,0  v!os,c,0,0  r  pM 

C 

2  fa,  0,0  rc,0,oYyios, i,0,0  ~vlos,c,0,o]^pn 

c 

2yT  n  /  \)  {  \) 

-^(vlos,s, 0,0  -vlos,c,0,0)Vpn) 

c 

7r,o0  —  CnYyios,s,0,0  —  '/  o,y ,  c ,  0 , 0  )^/j  ^  7,  / 
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or 
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Xv(i,n)=  Ar  exp 


2  co 


't, o  (  \  .  2/t, o  /  \2 

;  Vs, 0,0  rc, 0,0  j+  2  w>0,0  “  rc, 0,0  J 


— (fo-^fe.o.o-'i.o.ofe 

c 

7r,o(l  ,,s' ,o,o  ~  vios,c,o,o)TpnTsi 

c 

-  Cotrs,0,0  -  rcfifiYsi 


2®T,C 


(};los,s,  0,0  vlos,c,0,o}^ pn 


+ 

V  c 


i  (fi,0,0  G,0,0  'f\Vlos,s,(),0  vlos,c,0,o}^p 

c 

— 2^lv/oS>S,0,0  -vlos,c,0,0r\Tpnr 


(A51) 


We  observe  the  following. 

•  The  first  line  of  the  phase  expression  is  constant.  Therefore  it  does  not  contribute 
to  ascertaining  range  or  velocity  of  the  target. 

•  The  last  line  of  the  phase  expression  is  quadratic  in  index  n.  While  it  may 
defocus  the  Doppler  response,  it  does  not  contribute  to  its  location. 

•  The  second  and  third  lines  of  the  phase  expression  are  a  function  of  both  indices  n 
and  i.  They  represent  migration  terms.  They  are  typically  ignored  in  coarse 
resolution  GMTI  systems. 

Ignoring  these  terms  allows  us  to  simplify  the  video  phase  error  expression  to 

r  a  M 


Xv(i,n)=Ar  exp] 


2^IA(l  -  ^ofc.o.o  -  fi'Aokd 


c 


(yios,s,  0,0  vlos,c,0,o}^pft 


2  (fi',0,0  G ,  0 , 0  os ,  s  ,0,0  vlos,c,0,o}^p 


(A52) 


Typically  in  simple  GMTI,  the  first  phase  tenn  is  exploited  to  determine  target  range,  and 
the  second  phase  tenn  is  exploited  to  determine  target  Doppler,  although  the  third  term 
will  perturb  the  answer  somewhat.  The  third  term  is  a  result  of  the  RVPE  tenn 
previously  discussed. 
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We  identify  the  velocity  difference  as 


v/os,s,0,0  vlos,c,0fi > 


*5  (?o  +C,o,o) 

rc(t 

o  +  C,o,o ) 

v  r.v(?0  +C,o,o]| 

rc(t 

0  +tcfi,Q>\; 

vlos,t,  0,0 


(A53) 


where  all  values  are  instantaneous  values  at  time  corresponding  to  n  =  0  and  i  =  0 ,  that 
is  for  the  central  sample  of  the  central  pulse  of  the  CPI. 

Note  that  the  dot  product  is  merely  the  difference  between  line-of-sight  velocities  to  the 
target  location  and  the  SRP. 

Note  also  that  the  best  we  can  extract  from  the  data  is  the  quantity  {vtos  s_0  ()  -  vlos  c  0  0 ) . 

This  leaves  a  fundamental  ambiguity  in  being  unable  to  separate  the  target’s  own  velocity 
from  the  velocity  manifestation  due  to  target  location.  Even  if  the  target  were  stationary, 
there  would  be  a  velocity  component  due  to  its  offset  from  the  SRP.  It  is  precisely  this 
phenomenon  that  is  exploited  by  SAR. 
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“Do  not  worry  about  your  difficulties  in  Mathematics.  I  can  assure  you  mine  are  still 

greater.  ”  —  Albert  Einstein 
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